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Abstract. A framework for the investigation of disordered quantum systems in thermal 
equilibrium is proposed. The approach is based on a dynamical model — which consists of 
a combination of a double-bracket gradient flow and a uniform Brownian fluctuation — that 
'equilibrates' the Hamiltonian into a canonical distribution. The resulting equilibrium state is 
used to calculate quenched and annealed averages of quantum observables. 



1. Introduction 

In the conventional treatment of quantum statistical mechanics there is a natural division 
between (a) the system under study, which is treated quantum mechanically and whose states 
are subject to thermal fluctuations, and (b) the Hamiltonian of the system, which is treated 
essentially classically and is held fixed. For some quantum systems, however, the Hamiltonian 
itself may fluctuate for one reason or another. Questions that interest us in this connection, in 
particular, are: "How can a randomly fluctuating Hamiltonian approach its equilibrium state?" 
and "What is the form of the equilibrium distribution, and how do we calculate observable 
expectation values in equilibrium?" The former is a question of a dynamical nature, whereas 
the latter is a question of a static nature. The purpose of the present paper is to propose an 
approach to address these questions. Specifically, we shall derive a dynamical model having 
the property that a given initial Hamiltonian evolves randomly — but isospectrally — in such a 
way that the associated density function on the space of isospectral Hamiltonians approaches a 
steady state distribution given by the canonical ensemble. Furthermore, we apply the resulting 
equilibrium state to calculate thermal expectation values of other observables, leading to new 
physical predictions in some limiting (quenched and annealed) cases. 
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Figure 1. Spin in fluctuating magnetic 
field. The space of pure states of a spin- 
2 particle, in external magnetic field B, 
is the surface of the Bloch sphere. In 
statistical theory of quantum mechanics 
the state is represented by statistical 
distributions of the pure states, whereas 
the magnetic field B that specifies the 
Hamiltonian is held fixed. What happens 
if the direction of the field B is itself 
subject to a small fluctuation? 



2. Approach to equilibrium 

In classical statistical mechanics the notion of a gradient flow plays an important role in 
describing the approach to equilibrium: A system immersed in a heat bath naturally tends 
to release its energy into the environment and thus approach its minimum energy state, and 
this tendency is characterised by a Hamiltonian gradient flow. An equilibrium state is attained 
when this flow is on average counterbalanced by thermal noise due to a random interaction 
with the bath. Here the magnitude of the noise is determined by the temperature of the bath. 
Accordingly, the idea we are going to introduce here is a gradient flow equation on the space of 
Hamiltonians with the property that the eigenstates of an arbitrary initial Hamiltonian Hq at 
time t = tend toward alignment with those of a reference Hamiltonian, denoted by G. Thus, 
G plays the role of the 'fixed' Hamiltonian in conventional quantum statistical mechanics. The 
eigenstates of Ht thus evolve toward those of G under the flow. By introducing a suitable noise 
term, we are able to characterise the approach to an equilibrium distribution. 
The dynamical model for characterising approach to equilibrium is given by 

^ = -\[H t ,[H t ,G]] + [H t ,u t ], (1) 

where A € H+. Here {oJt} denotes a skew-symmetric matrix of independent white noise terms, 
and [H t ,u)t] is the Lie bracket of these with H t (see also pQ). Hence [H t ,cot] is symmetric and 
linear in both Ht and ujf The term [Ht, [Ht,G]] gives rise to the aforementioned gradient flow 
in the space of Hamiltonians. The Hermitian matrix G plays the role of the 'Hamiltonian of 
the Hamiltonians' in the sense that G determines the motion in the space of Hamiltonians. In 
particular, we can regard the linear function tr(HG) on the space of the totality of Hermitian 
matrices H as representing the 'energy' function defined on that space. 

There is an invariant measure (stationary solution) associated with the evolutionary equation 
JT]). This is given by the canonical density: 

p{H) <xexp(- \tr{HG)), (2) 



which can be used as a new basis for studying quantum statistical mechanics. In units h = 1 
the coupling A has dimension [Energy] , and has the interpretation of representing the inverse 
temperature for the 'Hamiltonian bath' (not to be confused with the thermal bath in which 
the system, and possibly also the apparatus determining the Hamiltonian, is immersed). The 
canonical density ([2]) can also be derived by entropy maximisation subject to the constraint that 
the energy tr(HG) has a definite expectation value. 



3. The gradient flow: double bracket equation 

Let us begin by examining properties of the gradient term in ([1]). Specifically, consider the 
following dynamical equation for Hermitian matrices: 

H H 

-^l = -X[Ht,[H t ,G}}. (3) 

Note that in terms of the Hermitian matrix X = \[H, G] the double-bracket evolution Q can 
be rewritten as 

d H 

-M=iX[H,X\, (4) 



which formally is just the Heisenberg equation of motion. However, owing to the //-dependence 
of X the evolution is nonunitary. The Hamiltonians Hq and G are both assumed nondegenerate. 



The flow induced by ([3j) satisfies the following properties: (i) the evolution is isospectral, i.e. the 
eigenvalues of Hq are preserved; and (ii) the evolution gives the 'alignment' \im.t-,oo[Ht, G] = 0. 

We remark that the double bracket flow was first introduced in the context of magnetism 
(Landau-Lifshitz equation) [2]. In its modern form it was introduced by Brockett [3] and has 
been successfully applied to many areas, such as optimal control, linear programming, sorting 
algorithms, and dissipative systems (see references cited in [3]). 

In the case of a 2 x 2 Hamiltonian the gradient flow equation ([3]) can be solved 
straightforwardly. In this case we express the Hamiltonian in terms of the Pauli matrices: 



H t = \u t \ + ^vcr-n t , 
where ri£ = (xt,yt,zt)- Similarly for the reference Hamiltonian G we write 



G=hvl + i ucr-g 



for a unit vector 
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A calculation shows that the solution to 



(5) 



(6) 



reads 



1 i uq — z/tanh(o;i — cq) i/sech(u;i — cq) e 1<?io 
2 ' v sech (u>t — cq ) e 1 ^ uq+u tanh (u>t — cq ) 



(7) 



where uj = Xufi, and cq = tanh 1 (cos^o) an d are initial values [4J. Furthermore, the 
eigenvalues of H t are time-independent, and we have 



lim Ht 

t— >oo 



n 







v 

U + V 



(8) 



Thus, the Hamiltonian is asymptotically diagonalised in the G-basis. Observe that trHt and 
detift are conserved quantities. Therefore, the flow induced by §3$) for fixed initial values uo and 
| no | is confined to a two-sphere C, which is isomorphic to the state space of a two-level system. 
It follows that in two dimensions we have the equivalence of the Schrddinger and Heisenberg 
pictures, even though the dynamical equation is not unitary. Since no and |no| are constant, we 
fix these and focus our attention on C parameterised by the dynamical coordinates (6t, 4>t)- 




Figure 2. Gradient flow with unitary motion on the 
sphere C. The vector field generated by the unitarily- 
modified gradient-flow @ is plotted. The first term in @ 
generates a rotation around the G-axis (which is chosen 
to be the z axis here), while the second term generates 
geodesic flows toward the south pole. The axis no of the 
initial Hamiltonian Hq spirals around the G-axis g and is 
asymptotically aligned with the latter (south pole in this 
example) . 



We remark that the dynamical equation ([3]) can be modified to include a unitary term: 

dH t 



dt 



-i[H t ,G]-X[H t ,[H t ,G]} 



(9) 



without greatly affecting its physical characteristics. In the 2x2 example, the only change 
occurs in the phase so that instead of <j>t = 4>o we have <j>t = 4>q + \it. 



4. Elements of stochastic differential geometry 

We now wish to introduce a Brownian term into the deterministic flow ([3]). Specifically, we 
consider a uniform Brownian field on the isospectral subspace of the space of Hermitian matrices 
(the sphere C in the 2x2 case). Before we proceed, however, it will be useful to recall how 
stochastic motions can be defined on a manifold. The basic process we consider is the Wiener 
process {Wt} defined on a filtered probability space (O,^ 7 , P). Here O is the sample space, T 
is a (7-field on Q, and P is the probability measure. The filtration of T determines the causal 
structure of (O,^ 7 , P). This is given by a parameterised family {Tt}§<t<oo of nested cr-subfields 
satisfying T s C F% C T for any s < t < oo. We say that {Wt} is a Wiener process if: (a) 
Wo = 0; and (b) {Wt} is Gaussian such that W t +h — Wt has mean zero and variance \h\ (see 
[5]). A process {at} is said to be adapted to the filtration {Tt} generated by {Wt} if its random 
value at time t is determined by the history of {Wt} up to that time. 

If {at} is ^-adapted, then the stochastic integral M t = J <7 s dW s exists, provided that {at} is 
almost surely square-integrable. If the variance of {Mt} exists, then {Mt} satisfies the martingale 
conditions E[|Mt|] < oo and E[A/t|jT s ] = M s , where E[— ] denotes expectation with respect to 
the measure P. The latter condition implies that given the history of the Wiener process up to 
time s the expectation of Mt for t > s is given by its value at s. 

A general Ito process is defined by an integral of the form 

xt = x + / /^ds + / a s dW s , (10) 
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where {fit} and {at} are called the drift and the volatility of {xt}- A convenient way of expressing 
(TlO]) is to write dxt = fj.tdt + atdWt, and to regard the initial condition xo as implicit. In the 
special case nt = ^{xt) and a t = a(xt), where /i(x) and a{x) are prescribed functions, the 
process xt is said to be a diffusion. 

This analysis can be generalised to the case of a diffusion {xt} taking values on a manifold 9JI, 
driven by an m-dimensional Wiener process {W t i }j = i v .. jm . Let V a be a torsion-free connection 
on dJt such that for any vector field £ a its covariant derivative in local coordinates is 

te(v b n = ^F + r^ c , (ii) 

where <5„ is the standard coordinate basis in a given coordinate patch. Suppose we have an Ito 
process taking values in 9Jt. Let xf denote the coordinates of the process in a particular patch. 
Then writing h ah = afa ht we define the drift process fi a by 

H a dt = dx a + iTg c /i bc dt - afdWl (12) 

Alternatively, we can write the covariant Ito differential as dx a = 5®(dx a + ^T^ c h hc dt), where 
is the dual coordinate basis. Then (fT2l) can be represented as 



dx a = n a dt + afdWi (13) 

If jJi a {x) and af(x) are m + 1 vector fields on SOT, then the general diffusion process on 971 is 
governed by a stochastic differential equation dx a = [i a (x)dt + o""(x)dW/, where dx a is the 
covariant Ito differential associated with the given connection (see Hughston [6]). 

For the characterisation of the diffusion process it suffices to specify a connection on 9JT, and 
a metric is not required. The quadratic relation dx a dx b = h ab dt, where h ab = afa bl , follows 
from the Ito identities dt 2 = 0, dtdWl = 0, and dW/ dW/ = d^dt. Then for any smooth function 
4>(x) on 971 we define the associated process (ft = 4>(xt), and Ito's formula takes the form 

dft = {V a ^)dx a + \{V a V b f)dx a dx b 

fi a V a f + \h ab V a V b A dt + afVafdWf (14) 



The probability law for xt is characterised by a density function p(x,t) on DJt that satisfies the 
Fokker-Planck equation 

^ = -V a (//» + \V a V b (h ab p). (15) 

The diffusion is said to be nondegenerate if h ab is of maximal rank. If g ab is a Riemannian metric 
on Wl and V a is the associated Levi-Civita connection, then if h ab = <J 2 g ab , the process x t is a 
Brownian motion with drift on 9Jt, with volatility parameter a. 



5. Diffusion model for thermalisation 

Consider a stochastic differential equation of the form 

dx a = p a dt + KafdW} (16) 

on a real manifold 9Jt. Here k is a constant, the drift p a is a vector field on 9Jt, and the vectors 
{erf} constitute an orthonormal basis in the tangent space of 9Jt. In this case the associated 
Fokker-Planck equation reads 

^ ptix) = -V a (/x» + ±/c 2 VV (17) 

For our model we require that the drift vector p a represent the double-bracket gradient flow (|3j) . 
This is achieved by choosing 

p a = — \ K 2 W a G, (18) 

where G{x) is a function on 971 given by tv(HG). Then it follows that there exists a unique 
stationary solution to (|17p . given by the canonical density 

J m exp{-XG{x))dV 

If 9Jt is the space of pure states, then we have a model for thermalisation of quantum states 
introduced by Brody & Hughston [7]. 



6. Two-dimensional case in more detail 

To illustrate these results in more explicit terms we consider a system consisting of a single 
spin-^ particle immersed in an external magnetic field. The Hamiltonian is then H = — B • S, 
where B denotes the field and S the spin vector. The direction of the field B, however, is subject 
to fluctuations around its stable direction, specified by G (directed along the z-axis). For the 
dynamical equation we obtain [1]: 



dd t = uo sin# t dt + v^dWj + dW t 2 ) 
d<t>t = -^ n r t V^(dW t 1 -dW?). 



(20) 



The associated Fokker-Planck equation reads 

p = -Lu(cose + sm6d e )p + 2v(d 2 e + d£)p, (21) 
where dg = d/d6 and = d/d(p. The asymptotic solution is the canonical density function: 

P( 9 >® = o -^U ^ exp (4Vcos#) . (22) 
27rsmh(^A / u) 



Direct substitution shows that (I22|) is the stationary solution to (121ft . 

It follows from ([22]) and the use of the volume element = \ sin 9d9d(j) that the equilibrium 
mean Hamiltonian is 

{ H) = l( U ° + ^ OS ^ f m )■ (23) 

\ tt - z/(cos0)a, y 

where 

^ ^ A XpL tanh(|A//) ^ 

We may regard the parameter A as representing the 'inverse temperature' for the Hamiltonian. 
If the noise level is high (A <C 1), then the direction of the external field B on the average lies 
close to the xy-plane so that {cos9)\ ~ 0. If the noise level is low A S> 1, then the field B on 
the average is parallel to the z-axis and we have {cos9)\ ~ — 1. We plot {cos9)\ as a function 
of t = 1/A (see Fig. ED . 



Figure 3. The plot of the expectation 
(cos9)\ as a function of the inverse 
Hamiltonian temperature r = 1/A. For 
r « 1 we have (cos9}\ ~ —1, whereas 
for r>lwe find (cos 9)\ ~ 0. 



7. Quantum statistical mechanics of disordered systems 

Now we consider how the statistical theory of Hamiltonians presented above can be applied to 
quantum statistical mechanics, when the system and the apparatus specifying the Hamiltonian 
are both immersed in a heat bath with inverse temperature (5. In this context it is natural to 
borrow ideas from the spin glass literature. We may take the averaged Hamiltonian {H)\ as the 
starting point of the analysis — this gives the analogue of an annealed average: 

tr(Oe-«">») , , 

{0) « = trVw.) (25) 

of an observable O. Such an averaging, however, will change the eigenvalues of H. 

Alternatively, we may use the 'unaveraged' Hamiltonian to compute the expectation of an 
observable O, and then take its average: 

/ tr(Oe~^) \ 

{0)Q = \ tr(e-^) A' (26) 

This gives the analogue of a quenched average. The canonical quenched average of the 
Hamiltonian G = a z is 



whereas the canonical annealed average of G is 



(G)a = A/^tanh 



tanh(| A/x) 



_2_ 

A/j 



(28) 



These results suggest a new line of studies on the extended quantum statistical mechanics of 
disordered systems. The plot below shows the annealed (blue) and quenched (red) averages of 
a z , as a function of the bath temperature T = 1/(3, for fixed A such that A~ x = 0.1. 
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Figure 4. Quenched and annealed av- 
erages of G. The functions (G)q and 
{G)a are plotted against the tempera- 
ture T = f3 , where we set A -1 = 0.1, 
v = 0, v = 1, and // = 2 so that G = 
a z . The 'quenched magnetisation' {<t z )q 
does not attain the maximum value 1.0 
at zero temperature unless A -1 = 0. 



8. Examination of higher-dimensional cases 

The geometry of higher-dimensional Hermitian matrices is somewhat more intricate than the 
two-dimensional case examined above. The space of N x N Hermitian matrices has the structure 
of the product x CP^ x CP^ -1 x • • • x CP 1 , where CP fc denotes the complex projective k- 
space. This space is considerably larger than the space CP^ of pure states upon which N x N 
Hermitian matrices act, and as a consequence the equivalence of the Schrodinger and Heisenberg 
pictures for a nonunitary motion would in general be lost. 

A 3 x 3 Hermitian matrix H can be written as H = E 1 \E 1 )(E 1 \ + E 2 \E 2 ){E2\ + E 3 \E 3 ){E 3 \, 
where {E{\ are the eigenvalues and {|i%)} are the associated eigenstates. The degrees of freedom 
for the energy eigenvalues correspond to the open space R 3 ; the remaining degrees of freedom 
corresponding to CP 2 x CP 1 are encoded in the specification of the energy eigenstates. Writing 
{|<7i)} for the eigenstates of G, we can express \E\) in the form: 

= sin itfcos l<p\gi) + sin ±i?sin \<pS\g 2 ) + cos ±$e iv \g 3 }, (29) 

which determines a point in CP 2 . There is a CP 1 worth degrees of freedom left for \E 2 ): 

\E 2 ) = ^cos \a sin \-Q cos \ip — sin \a sin ^cpe^j \g\) 

+ ^cos \a cos |#sin ^e 1 ^ + sin cos ^(pe^^J \g 2 ) + cos \a sin ^-de^lg^). (30) 

The specifications of \E\) and \E 2 ) leave no further freedom left for \E 3 ) and we have 

\E 3 ) = ^sin |a cos |$ cos |y> + cos \ct sin ^(pe 1 ^ \g\) 

+ ^sin \a cos \-d sin \<pe^ — cos \a cos ^tpe^^J \g 2 ) + sin \ a sin ^i9e 1T? | g 3 ). (31) 

In this manner we obtain the nine parameters required for the specification of an arbitrary 3x3 
Hermitian matrix. 



It should be remarked that the foregoing procedure is merely an example of how one might 
parameterise a generic N xJV Hermitian matrix in a systematic manner; the scheme is somewhat 
unconventional in that a generic 2x2 Hermitian matrix in this parameterisation is written 

/ 1 \ / cos<9 sinfle-^ \ ,„_,. 

H = "+{0 1 ) + "-{ S in8e^ -cose J' ^ 

where uj± = \{E\ ± E2). This set of coordinates is nevertheless convenient because it isolates 
invariant quantities u± from the coordinates (#, 4>) of the isospectral submanifolds. 

For our application in quantum statistical mechanics of disordered systems we are required 
to determine the partition function 

Z[X) = j e~ Xtl{HG) dV. (33) 

In the case of a 3 x 3 Hamiltonian, since the dynamical equation ([1]) preserves the three 
eigenvalues of H, the dynamical motion stays on the isospectral submanifold spanned by the six 
angular variables 1?, 99, a S [0, ir] and £, 77, f3 € [0, 2ir], with volume element 

dV = j^g sin a sin — cos $) sin cp da dj3 dt? dip d£ drj. (34) 

Since we are working in the G-basis, the calculation of tr(HG) is straightforward, and the integral 
(|33|) reduces to an expression analogous to the integral representation for a Bessel function. 

The choice of parameterisation adapted here for a generic N x N Hermitian matrix need not 
be the most adequate for our purpose. An alternative approach is to write H = U^EU, where 
E is a diagonal matrix with eigenvalues {Ei}, and U is a unitary matrix. Then the calculation 
of the partition function becomes semi-Gaussian, and we are left with an integration over the 
invariant Haar measure [8] . Ideas from matrix theory or random matrices might prove useful in 
the statistical analysis of disordered quantum systems introduced here. 
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